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Abstract. We report recent results from a multiphase transport (AMPT) model on 
the azimuthal anisotropics of particle momentum distributions in heavy ion collisions 
at the Relativistic Heavy Ion Collider. These include higher-order anisotropic flows 
and their scaling, the rapidity dependence of anisotropic flows, and the elliptic flow of 
charm quarks. 



PACS numbers: 25.75.Ld, 24.10.Lx 
1. Introduction 

Anisotropic flows in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) 
are sensitive to the properties of produced matter. This sensitivity not only exists 
in the larger elliptic flow Pfl El IH El El Ej but also in the smaller higher-order 
anisotropic flows [HI El El- Experimentally, scaling relations among hadron anisotropic 
flows have been observed jTTj, and according to the quark coalescence model they 
are related to similar scaling relations among parton anisotropic flows JU]- Also, 
anisotropic flows measured at flnite pseudorapidities are seen to depend strongly on 
the rapidity [HI [T21 UHl HH UHl HEj, and this has so far not been reproduced by 
theoretical models ^71 ^1 Q^l 120] ■ For charmed mesons, the observed transverse 
momentum spectra of electrons from their decays are found to be consistent with 
both limiting scenarios of perturbative QCD spectra without flnal-state interactions 
and complete thermalization including transverse expansion [21]. The charm elliptic 
flows given by the quark coalescence model show, however, marked differences between 
these two scenarios |22|. These interesting phenomena have recently been studied in 
a multiphase transport (AMPT) model that includes both initial partonic and flnal 
hadronic interactions [23 121I2S1 123 12ZI I2H1 • In this talk, we briefly review the AMPT 
model and discuss the results obtained from this model. 
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2. The AMPT model 

The AMPT model is a hybrid model that uses minijet partons from hard processes 
and strings from soft processes in the Heavy Ion Jet Interaction Generator (HIJING) 
model |22j as the initial conditions for modeling heavy ion collisions at ultra-relativistic 
energies. Time evolution of resulting minijet partons is then described by Zhang's parton 
cascade (ZPC) |3p| model. At present, this model includes only parton-parton elastic 
scatterings with an in-medium cross section derived from the lowest-order Born diagram 
with an effective gluon screening mass taken as a parameter for fixing the magnitude 
and angular distribution of parton scattering cross section. After minijet partons stop 
interacting, they are combined with their parent strings, as in the HIJING model with 
jet quenching, to fragment into hadrons using the Lund string fragmentation model as 
implemented in the PYTHIA program The final-state hadronic scatterings are 

then modeled by a relativistic transport (ART) model 

Since the initial energy density in Au + Au collisions at RHIC is much larger 
than the critical energy density at which the hadronic matter to quark-gluon plasma 
transition would occur |23 EH] , the AMPT model has been extended by converting the 
initial excited strings into partons [H]. In this string melting scenario, hadrons, that 
would have been produced from string fragmentation, are converted instead to valence 
quarks and/or antiquarks. Interactions among these partons are again described by the 
ZPC parton cascade model. Since there are no inelastic scatterings, only quarks and 
antiquarks from the melted strings are present in the partonic matter. The transition 
from the partonic matter to the hadronic matter is achieved using a simple coalescence 
model, which combines two nearest quark and antiquark into mesons and three nearest 
quarks or antiquarks into baryons or anti-baryons that are close to the invariant mass of 
these partons. The present coalescence model is thus somewhat different from the ones 
recently used extensively [Sll ESI 1301 IHI] for studying hadron production at intermediate 
transverse momenta. 

3. Anisotropic flows at midrapidity 

Using the AMPT model with string melting, we have studied anisotropic fiows in the 
pseudorapidity range |?7| < 1.2 in minimum bias Au + Au collisions at ^/s = 200 AGeV. 
In Fig. ^ we show the final anisotropic fiows f„ of charged hadrons, defined by the 
average < cos(n0) > of the azimuthal distributions of their transverse momenta, as 
functions of transverse momentum pt for parton scattering cross sections ap = 3 and 
10 mb, together with recent experimental data from the STAR collaboration ^11^. It 
is seen that the parton scattering cross section ap = 3 mb underestimates the data at 
higher px (> 1 GeV/c) while cXp = 10 mb seems to give a better fit to the data. The 
values of hadronic vq are in agreement with the data within error bars, although they 
are essentially zero. The V4 of charged hadrons exhibits a stronger sensitivity to the 
parton cross section than their V2, especially at higher pt, and is thus a more sensitive 
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Figure 1: Anisotropic flows V2 (a), ^4 (b), 
and Vq (c) of cliarged liadrons in \ri\ < 1.2 
from minimum bias Au + Au collisions at 
y/s = 200 AGeV as functions of px for 
parton cross sections of 3 (open squares) 
and 10 (solid squares) mb. Data are from 
the STAR Collaboration 11 . 
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Figure 2: Transverse momen- 

tum dependence of midrapidity parton 
anisotropic flows V2, V4 and Vq from 
minimum bias events for Au + Au at 
y/s = 200 AGeV with parton scatter- 
ing cross section 10 mb. Also plotted 
are t;| (open triangles) and V2 (open 
stars). 



probe to the initial partonic dynamics in relativistic heavy ion collisions. 

The experimental data shown in Fig. ^ on higher-order anisotropic flows show 
the scaling relation Vnipr) ~ '^2^'^{Pt) [H]- In the naive quark coalescence model |SZj 
that only allows quarks with equal momentum to form a hadron, meson higher-order 
anisotropic flows v^^m and vq^m are related to parton anisotropic flows Vn,q by |21 El IHI] 

Vi,M{PT) 1 1 ^^4,9 (Pt/2) V&,m{Pt) 1 ( V4,q{pT/2) V6^q{pT/2) \ 

vIm{Pt) 4 ^2vUpT/2y vIUpt) 4 [vlipT/2) + vl^ipT/2)) ' ^ > 
The meson higher-order anisotropic flows thus satisfy scaling relations if such relations 
exist among quark higher-order anisotropic flows. 

In Fig. 121 we show the pt dependence of the anisotropic flows V2, and of 
midrapidity partons obtained from the AMPT model with string melting and a parton 
scattering cross section of 10 mb for above reaction. Also shown in Fig. |21are f| (open 
triangles) and (open stars). Comparing them to ^4 and fg, respectively, indeed show 
that the parton anisotropic flows satisfy the scaling relation Vn^qipr) ~ '^2^q{PT)- With 
a parton scaling factor of 1, the naive coalescence model would lead to the following 
scaling relations for meson anisotropic flows: 

^4,M(pr) _ 3 vq,m{Pt) _ 1 
vImM^ vIm{Pt)^2- 

The resulting hadron scaling factors of 3/4 and 1/2 are, however, smaller than those 
(~ 1.2) extracted from measured anisotropic flows of charged hadrons. Since the naive 
quark coalescence model does not allow hadron formation from quarks with different 
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momenta as in more realistic quark coalescence models jSH ESI ISS]; it is not expected 
to give a quantitative description of the experimental observation. Such effects are, 
nevertheless, included in the AMPT model, which have been shown in Fig. ^ to 
reproduce measured hadron anisotropic flows. 

4. Pseudorapidity dependence of anisotropic flows 
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Figure 3: Pseudorapidity dependence 
of Vi from minimum bias events of Au 
+ Au collisions at = 200 AGeV in 
the string melting scenario with parton 
scattering cross sections dp = 3 (open 
squares) and 10 (solid squares) mb 
as well as the scenario without string 
melting (solid triangles). Data are from 
the STAR collaboration [TT]- 



Figure 4: Pseudorapidity dependence 
of V2 from minimum bias Au + Au 
colhsions at ^/s = 200 AGeV in the 
string melting scenario with parton cross 
sections of 3 (open squares) and 10 (solid 
squares) mb and the scenario without 
string melting (solid triangles). Data are 
from the PHOBOS (stars) [13] and STAR 
collaborations (solid circles) jl5j . 



Results from the AMPT model on the pseudorapidity dependence of vi for charged 
hadrons from minimum bias events of Au + Au collisions at ^/s = 200 AGeV are 
shown in Fig. 01 for the scenarios of string melting with parton scattering cross sections 
CTp = 3 (open squares) and 10 mb (solid squares) as well as for the scenario without 
string melting (solid triangles). Also included in Fig. El are recent data from the STAR 
collaboration (solid circles) [TT]. Both scenarios can reproduce approximately the data 
around the mid-pseudorapidity region, i.e., Vi is flat (essentially zero) around mid-?]. 
For vi at large {rjl, the string melting scenario with both parton scattering cross sections 
CTp = 3 mb and 10 mb underestimates significantly the data. On the other hand, the 
scenario without string melting seems to give a good description of Vi at large \7]\. 

The predicted pseudorapidity dependence of charged hadron V2 from the same 
reaction is shown in Fig. 0] together with preliminary data from the PHOBOS 
collaboration (solid stars) [13] and the STAR collaboration (solid circles) The 
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string melting scenario with cXp = 10 mb (solid squared) is seen to describe very well the 
data for V2 around mid-rj {\r]\ < 1.5) but overestimates the data at large pseudorapidity. 
The overestimation of may stem from the use of constant parton scattering cross section 
in the AMPT model. Since the properties of partonic matter at different rapidities may 
not be the same in heavy ion collisions at RHIC, different parton cross sections may 
have to be used. Comparison between theoretical results and the experimental data for 
elliptic flow indicates that a larger dp = 10 mb is needed at midrapidity but a smaller 
(Tp = 3 mb (open squares) gives a better description at large pseudorapidity. Also shown 
in Fig. lUare results obtained from the scenario without string melting (solid triangles), 
and they are seen to also describe the data at large pseudorapidity (|?7| > 3). Therefore, 
the scenario without string melting can describe simultaneously the data on vi and f 2 
at large pseudorapidity {\ri\ > 3). These interesting features imply that initially the 
matter produced at large pseudorapidity {\ri\ > 3) is dominated by strings while that 
produced around mid-rapidity {\r]\ < 3) mainly consists of partons. This is a reasonable 
picture as particles at large rapidity are produced later in time when the volume of the 
system is large and the energy density is small. 



5. Charm flows at RHIC 




Figure 5: Transverse momentum dependence of V2 for different parton flavors in 
midrapidity from minimum bias Au + Au collisions at ^/s = 200 AGeV. Results are 
for parton cross sections fip = 3 (a) and 10 mb (b). 



Using the AMPT model in the string melting scenario with same parton scattering 
cross section of dp = 3 or 10 mb for all quarks, we have also studied the pr-dependence 
of V2 for partons of different flavors in the midrapidity from minimum bias Au + Au 
collisions at ^/s = 200 AGeV. Using the current quark mass of 10 MeV for d quark, 6 
MeV for u quark, 200 MeV for s quark, and 1.35 GeV for c quark, the results are shown 
in Fig. El R is seen that the quark V2 increases with increasing pt, and at pt < 1-5 
GeV/c it is smaller for the heavier c quark than for the lighter d and s quarks. This mass 
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dependence of f 2 at low px is similar to the mass ordering of hadron elliptic flows in the 
hydrodynamic model which assumes that the matter is in local thermal equilibrium and 
thus develops a large collective radial flow. However, instead of continuing increase of V2 
with respect to pt as in the hydrodynamic model, the parton V2 in the transport model 
reaches a maximum value at certain large pt, indicating that high momentum particles 
do not achieve thermal equilibrium with the bulk matter. These features are seen for 
both values of parton scattering cross sections, except that the value of f 2 is larger for a 
larger parton scattering cross section. As in the case of light hadrons, study of charmed 
hadron flows offers the possibility to understand the charm quark interactions in the 
quark-gluon plasma. 

6. Summary and discussions 

In summary, using the AMPT model, we have studied the partonic effects on anisotropic 
flows in heavy ion collisions at RHIC. We find that measured V2, f4 and Vq of charged 
hadrons at midrapidity in Au + Au collisions at y/s = 200 AGeV can be described by a 
parton scattering cross section of about 10 mb and that V4 is a more sensitive probe to 
the initial partonic dynamics in these collisions than V2- Moreover, higher-order parton 
anisotropic flows are nonnegligible and satisfy the scaling relation Vn^qipr) ~ "^2,9^ (Pt), 
which leads naturally to the observed similar scaling relation among hadron anisotropic 
flows if the partonic matter hadronizes via the coalescence of quarks and antiquarks. The 
results on the rapidity dependence of anisotropic flows suggest that a partonic matter 
is formed during early stage of relativistic heavy ion collisions only around midrapidity 
and that strings remain dominant at large rapidities. We have also studied the elliptic 
flow of charm quarks by using different scattering cross sections, and find that it has a 
maximum value similar to that for light quarks and is also sensitive to the charm quark 
scattering cross sections. The charm quark elliptic flow shows, on the other hand, a 
very different dependence on the transverse momentum due to its large mass. 

In reproducing the experimental data on hadron anisotropic flows, a large parton 
cross section of about 10 mb is needed in the AMPT model. Comparing to the cross 
section given by the perturbative QCD, this value is an order of magnitude larger, 
indicating that nonperterbative effects are important in the quark-gluon plasma. Indeed, 
recent lattice gauge studies have shown that hadron spectral functions including those 
consisting of heavy charm quarks survive in the quark-gluon plasma SSj 39; . Including 
the effect of these quasi bound states is expected to enhance the parton-parton scattering 
cross section. 
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